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ABSTRACT 

We report the discovery of the second and third pulsating extremely low mass white dwarfs ( WDs) , 
SDSS J111215.82+111745.0 (hereafter J1112) and SDSS J151826. 68+065813.2 (hereafter J1518). 
Both have masses < 0.25 Mq and effective temperatures below 10, 000 K, establishing these pu- 
tatively He-core WDs as a cooler class of pulsating hydrogen- atmosphere WDs (DAVs, or ZZ Ceti 
stars). The short-period pulsations evidenced in the light curve of J1112 may also represent the first 
observation of acoustic (p-mode) pulsations in any WD, which provide an exciting opportunity to 
probe this WD in a complimentary way compared to the long-period (/-modes also present. J1112 is a 
T cS — 9590 ±140 K and log g = 6.36 ±0.06 WD. The star displays sinusoidal variability at five distinct 
periodicities between 1792 — 2855 s. In this star we also see short-period variability, strongest at 134.3 
s, well short of expected g-modes for such a low-mass WD. The other new pulsating WD, J1518, is a 
T cfr = 9900±140 K and log 5 = 6.80±0.05 WD. The light curve of J1518 is highly non-sinusoidal, with 
at least seven significant periods between 1335 — 3848 s. Consistent with the expectation that ELM 
WDs must be formed in binaries, these two new pulsating He-core WDs, in addition to the prototype 
SDSS J184037. 78+642312. 3, have close companions. However, the observed variability is inconsistent 
with tidally induced pulsations and is so far best explained by the same hydrogen partial-ionization 
driving mechanism at work in classic C/O-core ZZ Ceti stars. 

Subject headings: binaries: close - Galaxy: stellar content - - Stars: individual: SDSS 
J111215.82+111745.0, SDSS J151826.68+065813.2 — Stars: white dwarfs — vari- 
ables: general 



1. INTRODUCTION 

There are many pulsation instability strips on the 
Hertzsprung-Russell diagram, including the DAV (or ZZ 
Ceti) instability strip, driven by a hydrogen partial ion- 
ization zone in the cool, hydrogen atmosphere (DA) 
white dwarfs (WDs). Seismology using the non-radial 
gravity-mode (g-mode) pulsations of DAVs has enabled 
us to constrain the mass, core and envelope composi- 
tion, rotation rate, and the behavior of convectio n in 
these objects (see reviews by IWinget fc Ke^l2008l and 
iFontaine fc Brassard 200§|). However, the roughly 150 
DAVs known to date have masses between 0.5 — 1.1 Mq, 
implying that they all likely contain C/O-cores. Low- 
mass, He-core WDs are likely to pulsate as well, in- 
cluding extremely low -mass (ELM, < 0.25 Mq) WDs 
(jSteinfadt et alJl20Toh . 

ELM WDs are the byproducts of binary evolution. 
Mass loss during at least one common-envelope phase 
has precluded the ELM WD from igniting He in its core, 
and the result is an underweight WD with a core de- 
void of C/O. These WDs have been known for some 
time as companions to pulsars, but recently their num- 
bers have grown dramatically as a result of the ELM 
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Survey, a targeted spectroscopic search for ELM WDs 
(IBrown et al.ll20Tol: iKilic et aL|[201lt iBrown et alj|20ll 
IKilic et al.H2012D . 

Leveraging the numerous new low-mass WDs found 
by the ELM Survey, the first putatively He-core ELM 
WD (SDSS J184037.38+6423 12.3, hereafter J184 0) was 
recently discovered to pulsate (|Hermes eTaT]|20"l^ ). This 
~0.17 Mq WD varies at a dominant period of roughly 
4698 s, with a high-amplitude (> 5%), non-sinusoidal 
pulse shape. The discovery of pulsating ELM WDs pro- 
vides the first chance to apply the tools of asteroseis- 
mology to these low-mass, presumed He-core WDs; we 
have an exciting opportunity to prob e the internal struc- 
ture of these exotic comp act objects (jCorsico et al.ll2012t 
IVan GrooteUtalllMl . 

Here we announce the discovery of pulsations in two 
new ELM WDs: SDSS J111215.82+111745.0 (here- 
after J1112) and SDSS J151826.68+065813.2 (hereafter 
J1518). These objects bring to three the number of pul- 
sating ELM WDs known, further populating the empiri- 
cal instability strip of ELM WDs. These objects belong 
to an extension of the instability strip from the C / O-core 
DAVs, and are most likely driven by the same mechanism 
acting in the classical ZZ Ceti stars. 

In Section [5] we detail our discovery of pulsations in 
J1112. We outline our new spectroscopic and photo- 
metric observations of this ELM WD, and provide a ro- 
bust identification of its highest-amplitude periodicities. 
We also detail evidence for short-period variability which 
may be evidence of acoustic (p-mode) pulsations in this 
ELM WD. In Section [3] we describe the discovery of vari- 
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ability in the ELM WD J1518 and list the significant 
periods determined from our relatively sparse coverage 
so far. We conclude with a discussion of these discover- 
ies in Section 01 

2. SDSS J1112+1117 

2.1. Spectroscopic Observations 

iBrown et al.1 ()2012h presented a preliminary fit of T e g 
= 9400 ± 490 K and logg = 5.81 ± 0.12 based on a 
single spectrum of this go — 16.2 ma^3 WD from the 
FLWO 1.5 m telesco pe using the FAST spectrograph 
(|Fabricant et al.lll998f) . We have obtained an additional 
26 spectra using the FLWO 1.5 m telescope and an ad- 
ditional 6 spectra us ing the Blue Channel Spectrograph 
(jSchmidt et al.lll989| ) on the 6.5m MMT. The time-series 
spectroscopy reveals that J1112 is short-period binary, as 
described below. 

2.1.1. Atmospheric Parameters 

We have phased (see Section 12. 1.2j) and co-added our 
new spectroscopic observations to determine the atmo- 
spheric parameters of the primary ELM WD visible 
in J1112. Our observations cover a wavelength range 
from 3700 — 4500 A. The model atmo spheres used for 
this a nalysis are described at length in iGianninas et all 
(|2011) and employ the new Star k broadening profiles 
from iTremblav fcB crgcron (2009). Models where con- 
vective energy transport becomes important are com- 
puted using the M L2/q = 0.8 prescripti on of the mixing- 
length theory fsee lTremblay et"al~ll2010( ). However, since 
we are dealing with ELM WDs, the model grid needed 
to be extended down to lower surface gravities. Thus the 
model grid used in this analysis covers the ranges in T g 
from 4000 to 30,000 K in steps ranging from 250 to 5000 
K and log g from 5.0 to 8.0 in steps of 0.25 dex. 

The method used for fitting the observations relies 
on the so-called spect roscopic technique, described in 
IGianninas et al.l f|201 lh and references therein. The main 
difference here is that we fit the Balmer lines we observe 
up to and including H12, as the considerably lower sur- 
face gravity means that these higher Balmer lines are 
still observed. Since the higher Balmer lines are sensi- 
tive m ostly to log g (see Figure 2 of lTremblav fc Bergeron! 
2009) , the inclusion of these extra lines further constrains 
our measurement of the surface gravity 

The uncertainty for each parameter is calculated by 
combining the internal error, which is the dominant 
source of uncertainty, obtained from the covariance ma- 
trix of the fitting algorithm with the external error, ob- 
tained from multiple observations of the same object, 
estimated for DA stars at 1.2 % in T e g and 0.038 dex in 
log g (see iLiebert et al.l 120051 for details) . The effect on 
the uncertainties caused by different values of S/N are 
included in the internal error. 

Our final fit to the phased and co-added spectra of 
J1112 is shown in the top panel of Figure [1] and yield 
T cff = 9590 ± 140 K and log .9 = 6.36 ± 0.06. This corre- 
sp onds to a mass of ~ 0.17 Mq using the He-core models 
of lPaneiet all ((20071) . 

In addition to the Balmer series, the Ca II K line at 
3933 A is also observed in the spectra of J1112. For the 
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Fig. 1. — Spectroscopic observations of J1112. The top panel 
shows the summed and phased spectra, with a model fit to the 
H7-HI2 lines of the Balmer series. This model derives the primary 
parameters in Section [2.1.1l The middle panel shows our new radial 
velocity observations of J1112 over six epochs, and the bottom 
panel shows those data phased to the orbital period of 4.13952 hr. 

purposes of this analysis, we simply exclude the wave- 
length range where that metal line is present so that it 
does not affect either the normalization of the individ- 
ual Balmer lines nor the actual fits themselves. This Ca 
feature phases with the radial velocity variations of the 
Balmer lines, and is thus not interstellar but rather in 
the photosphere of the WD. We will not discuss this Ca 
feature further in this paper. 

2.1.2. Radial Velocity Observations 

The bottom two panels of Figure Q] show the radial ve- 
locity curve of the 32 spectra we have obtained of J1112. 
We compute the orbita l elements of this s i ngle-l ined bi- 
nary using the code of iKenvon fc Garcial (|1986l ). which 
weights each velocity measurements by its associated er- 
ror. Our spectroscopic observations find that J1112 is 
in a 4.13952 ± 0.00024 hr orbital period binary with a 
K = 116.2 ± 2.8 km s _1 radial velocity semi-amplitude. 

This yields a mass function of / = 0.028 ± 0.003 Mq, 
which constrains the minimum mass of the unseen com- 
panion to M2 > 0.14 Mq assuming a 0.17 Mq primary. 
The nature of this companion has no direct bearing on 
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Fig. 2. — The top panel shows high-speed photometry of J1112 
from a representative run on 2012 February 1. The brightest com- 
parison star is shown in blue, offset by —8%. Average point-by- 
point photometric errors are also shown. The bottom panel shows 
a Fourier transform of our entire data set to date, some 70.2 hr 
of observations from 2012 January to 2012 April. We also display 
in red the Fourier transform of the residuals after prewhitening by 
the five highest-amplitude periods listed in Table [2] and mark those 
periods with green tick marks at the top of the panel. We mark 
the 4(A) significance level as a dashed green line. 

the discovery of pulsations in the primary ELM WD, but 
it may be contribute enough flux to partially contaminate 
the spectral fits of the primary. It also suggests that un- 
less the system is inclined to less than 30 degrees (which 
would happen at random less than 15% of the time), 
the companion is most likely another He-core WD with 
a mass below 0.45 M Q . Similar single- lined ELM WDs 
with likely He-core WD companions have be en found in 
the E LM Survey (i.e., SDSS J1005+3550 in lKilic et all 

i2on . 

We have folded our photometric observations on this 
orbital period and do not see evidence of eclipses, at a 
limit of 0.5%. This suggests that the inclination of the 
system is i < 80 degrees, although this does not signifi- 
cantly constrain the nature of the unseen companion. 

2.2. Photometric Observations 

We obtained high-speed photometric observations of 
J1112 at the McDonald Observatory over four months, 
from 2012 January to 2012 April, for a total of more 
than 70.2 hr of coverage. We used the Argos instrument, 
a frame-transfer CCD mount ed at the prime focus of th e 
2.1m Otto Struve telescope (jNather fe Mukadaml [20040 . 
to obtain 5 — 10 s exposures on J1112; a full journal of 
observations can be found in Table [TJ Observations were 
obtained through a 2mm BG40 filter to reduce sky noise. 

We performed weighted, circular, aperture pho- 
tometry on the calibrated frames using the external 
IRAF package ccd hsp written by Antonio Kanaan 
ijKanaan et al.l 12002ft . We divided the sky-subtracted 



TABLE 1 

Journal of photometric observations. 



Run 


UT Date 


Length Seeing 


Exp. 








(hr) 


(") 


to 






SDSS J1112-I-1117 






A2568 


2012 


Jan 21 


1.9 


3.3 


5 


A2571 


2012 


Jan 22 


1.3 


2.2 


10 


A2574 


2012 


Jan 23 


0.8 


2.0 


5 


A2577 


2012 


Jan 24 


2.2 


1.3 


10 


A2578 


2012 


Jan 25 


2.5 


3.4 


10 


A2581 


2012 


Jan 27 


6.5 


3.0 


5 


A2586 


2012 


Jan 29 


2.6 


1.7 


5 


A2588 


2012 


Jan 30 


3.2 


2.0 


5 


A2590 


2012 


Jan 31 


3.0 


1.2 


5 


A2594 


2012 


Feb 1 


3.5 


2.1 


5 


A2597 


2012 


Feb 2 


2.6 


1.9 


10 


A2600 


2012 


Feb 3 


2.8 


1.5 


5 


A2605 


2012 


Feb 15 


3.9 


1.6 


10 


A2607 


2012 


Feb 16 


2.4 


1.8 


10 


A2611 


2012 


Feb 19 


5.0 


2.1 


10 


A2614 


2012 


Feb 20 


2.9 


1.8 


10 


A2616 


2012 


Mar 12 


2.4 


2.4 


5 


A2618 


2012 


Mar 13 


2.9 


1.1 


10 


A2621 


2012 


Mar 14 


3.6 


1.1 


10 


A2638 


2012 


Mar 17 


3.1 


2.5 


10 


A2641 


2012 


Apr 14 


1.9 


3.3 


5 


A2656 


2012 


Apr 19 


3.7 


2.0 


5 


A2659 


2012 


Apr 20 


2.2 


1.5 


5 


A2662 


2012 


Apr 21 


1.7 


1.3 


5 


A2675 


2012 


Apr 24 


1.7 


1.1 


5 






SDSS J1518+0658 






A2626 


2012 


Mar 15 


0.7 


3.3 


5 


A2634 


2012 


Mar 16 


2.2 


2.2 


10 


A2639 


2012 


Mar 17 


3.4 


3.4 


10 


A2647 


2012 


Apr 16 


4.9 


1.7 


5 


A2650 


2012 


Apr 17 


1.7 


2.0 


5 


A2660 


2012 


Apr 20 


1.7 


1.2 


5 


A2664 


2012 


Apr 21 


3.9 


1.9 


10 


A2673 


2012 


Apr 23 


0.7 


1.5 


5 


A2676 


2012 


Apr 24 


0.4 


1.6 


10 


A2684 


2012 


Jun 18 


2.7 


2.1 


10 


A2686 


2012 


Jun 19 


2.9 


1.8 


10 


A2688 


2012 


Jun 20 


2.3 


2.4 


5 


A2705 


2012 


Jul 13 


3.1 


1.3 


5 



light curves by the brightest comparison star in the field, 
SDSS J111211. 51+111648.4 (g = 15.0 mag), to remove 
transparency variations, and applied a timing correction 
to each observation to account for the motio n of the 
Earth around the barycenter of the solar system (jStumpfll 
1980: iThompson fc MullalTvl[2009l ). 

The top panel of Figure [2] shows a portion of a typi- 
cal light curve for J1112, obtained on 2012 February 1, 
and includes the brightest comparison star in the field 
over the same period (the scatter is large because the 
other comparison stars in the field used to construct this 
divided light curve are all fainter than g — 16.9 mag). 
The bottom panel of this figure shows a Fourier trans- 
form (FT) utilizing all 38, 863 light curve points collected 
thus far. We display the 4(A) significance line at the bot- 
tom of Figure^ calculated from the average amplitude, 
(A), of an FT within a 1000 /iHz region in steps of 200 
/iHz, after pre-whitening by the five highest-amplitude 
periodicities. 

The pulse shape of J1112 appears quite sinusoidal, 
and is nearly solved with five independent periodicities. 
Those periods have been identified in decreasing order of 
amplitude in Table [5] For more realistic estimates, the 
cited errors are not formal least-squares errors to the data 
but rather the product of 10 5 Monte Carlo simulations 
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TABLE 2 

Frequency solution for SDSS J111215. 82+111745.0 



ID 


Period 


hi' ccjiic i ic v 


Amplitude 


S/N 




00 


(jiHz) 


(mma) 




h 


2258.528 ± 0.003 


442.7662 ±0.0007 


7.49 ± 0.08 


26.5 


h 


2539.695 ± 0.005 


393.7480 ±0.0007 


6.77 ± 0.09 


23.0 


h 


1884.599 ± 0.004 


530.6170 ±0.0011 


4.73 ± 0.08 


16.9 


U 


2855.728 ±0.010 


350.1734 ±0.0013 


3.63 ± 0.09 


11.5 


h 


1792.905 ± 0.005 


557.7542 ± 0.0017 


3.31 ± 0.08 


11.9 


h 


134.275 ± 0.001 


7447.388 ± 0.010 


0.44 ± 0.08 


U 


h 


107.56 ±0.04 


9297.4 ±3.6 


0.38 ± 0.14 


4.1 



Note. — 1 mma — 0.1% relative amplitude 

of perturbed data us ing the software package Period04 
(|Lenz fc Bregeiil2005h . The signal-to -noise calculation is 
based on the amplitude of the variability as compared 
to the average amplitude of a 1000 /iHz box centered 
around that variability, after pre-whitening by the five 
highest-amplitude periodicities. 

A simultaneous linear least-squares fit, fixing these five 
periods, shows that this variability is quite stable in both 
amplitude and phase. The /i periodicity is especially sta- 
ble in phase, with an r.m.s. scatter less than 7 s between 
our four months of data — better than 3% of the 2258.5 
s period. The phase stability of J1112 is reminiscen t 
of hot DAVs such as G117-B15A (|Kepler et all [20051) . 
and could be monitored long-term for periodic deviations 
in arriva l times to constrain any possible circumbinary 
planets (jMullallv et al.ll2008D . There is slightly more 
scatter about the amplitudes measured from month to 
month, which are more sensitive to variations in photo- 
metric conditions. The only periodicity with a consis- 
tently decreasing amplitude is /2, which showed an am- 
plitude of 0.7184% ± 0.0098% in 2012 January decrease 
to 0.573% ± 0.016% in 2012 April. 

None of these periods are an integer harmonic of the 
4.13952 hr (14902.27 ± 0.86 s) orbital period. Given our 
cited uncertainties, fs, the closest, is more than 8-a from 
8 x fori,. Thus, tidally induced pulsations cannot properly 
explain the observed multi-periodic variability. Instead, 
we conclude that these are global, non-radial <?-mode pul- 
sations driven to observa bility by the sam e mechanism at 
work in classical DAVs (|Brickhilllll991f) . The timescale 
of this variability is considerably longer than for the pul- 
sations seen in C/O-core DAVs. However, it is consistent 
with the expectation that the periods of pulsation modes 
roughly scale with the dynamical timescale for the whole 
star, n oc p -1 / 2 , and are thus much longer for the low- 
surface-gravity ELM WDs. 

2.3. Potential p-mode Pulsations 

In addition to the relatively high-amplitude, long- 
period variability observed in J1112, we see evidence 
for low-amplitude variability on much shorter timescales. 
These periodicities are included at the bottom of Tabled 
in decreasing order of S/N. This S/N value is conser- 
vative: We have not pre- whitened by the variability in 
question for this calculation, which effectively considers 
some signal as noise in this estimate. We have identified 
all periodicities with S/N > 4.0, and italicize these S/N 
values to indicate they were calculated in a different way. 

We have also computed the probability that each of 
the short-period detections is real by computing the 
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Fig. 3. — A zoom of the high-frequency regions in the FT of 
the entire J1112 data set that display evidence for short-period 
variability, potential p-mode pulsations at 134.3 s (top) and 107.6 
s (bottom). The original FT is shown in black, and the red shows 
the residuals after pre-whitening by the highest-amplitude peak. 
The dashed blue and green lines show the 3(A) and 4(A) lines, 
respectively. The lower panel shows the spectral window in blue 
centered around each periodicity. 

Fa lse Alarm Pro bability, using the formalism described 
in IKeplerl (|l99l . We find that f e and f 7 have a FAP 
> 99.9%. (There is an additional peak at 119.552 s with 
S/N = 3.7 and a FAP of 99.8%. However, this period- 
icity is sufficiently close to the 119.667 s periodic drive 
error of the 2.1 m Otto Struve telescope that we will not 
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include it in our formal frequency solution.) Figure |3] 
shows a zoomed-in portion of the FT around 134.3 s and 
107.6 s using our entire data set. 

This variability is coherent enough to reach significant 
amplitude over four months of observations. Some of our 
longest individual runs also evidence this variability, such 
as the 6.5 hr run on 2012 Jan 27 (with a 1.24 ±0.36 mma 
signal at 134.2 ± 2.2 s) and the 3.5 hr run on 2012 Feb 1 
(with a 1.17 ± 0.40 mma signal at 107.9 ± 1.8 s). These 
two peaks are also fairly significant if we use just our 2012 
January data: The 134.275 s mode has 0.55 ± 0.17 mma 
amplitude (> 99.9% FAP) while the 107.557 s mode has 
0.46 ± 0.14 mma amplitude (98.2% FAP). 

This variability is too short to be explained as g- 
mode pulsations without invoking implausibly high val- 
ues of the spherical harmonic degree. A recent non- 
ad iabatic pulsation an alysis relevant to low-mass WDs 
bv lCorsico et al.l (|2012| ) found that g- modes of low radial 
order (and thus the shortest period) are stable to pulsa- 
tions and should not be driven to observability. Their 
calculations found that unstable g- modes in a 0.17 M 
WD have radial orders k > 9 and periods longer than 
1100 s. Even if we ignore their conclusion that an £ = 1, 
k = 1 mode is stable, they find this low est radial order 
mode has a period ~249.5 s. Similarly, iSteinfadt et al.l 
( 2010) also found that an £ = 1, k = 1 .g-mode for a 0.17 
M Q WD has a ^245 s period. 

The short-period variability seen in J1112 is also incon- 
sistent with nonlinear combination frequencies present in 
the non-sinusoi dal light curves of many classical DAVs 
(Brickhill 1992). For one, the light curve of J1112 is ex- 
tremely sinusoidal. Additionally, the short-period vari- 
ability is not a multiple of any of the five low-frequency 
modes, nor is it a combination of different modes. 

Instead, we propose that this variability is caused by 
acoustic or pressure (jo-mode) pul s ations driven to ob- 
servability in J1112. iCorsico et al.l (|2012l ) find that low- 
order p-modes are pulsationally unstable, and have pe- 
riods ranging from 109 — 7.5 s for their 1 < k < 29 
models of a 0.17 M© He-core WD. The 134.3 s period 
we observe in J1112 is slightly longer than this predicted 
range, which suggests some uncertainty in identifying the 
true nature of these instabilities. Still, should these hold 
up as acoustic modes, this would mark the first detection 
of p-mode pulsations in any WD. We discuss the impact 
of this discovery in Section |4] 

3. SDSS J1518+0658 
3.1. Spectroscopic Observations 

iBrown et alJ (|2012h present the spectroscopic discovery 
data for this go = 17.5 mag WD from the Blue Channel 
spectrograph on the 6.5m MMT. They use 41 separate 
spectra over more than a year to determine the system 
parameters, and find that J1518 is in a 14.624 ±0.001 hr 
(52646.4±3.6 s) orbital period binary with a K = 172±2 
km s _1 radial velocity semi-amplitude. 

We have fit their 41 phased and co-added spec- 
tra with the extende d stell ar atmosphere models of 
iTremblay fc Bergeron] (|2009f ). as described in Sec- 
tion |2"TH This fit formally yields T cS = 9900 ± 140 
K and log g = 6.80 ± 0.0 5 for J1518, which corresponds 
to a mass of ~0.23 M© (|Panei et al.l 120071 ) . Given the 
mass function (/ = 0.322 ± 0.005 M©), the minimum 
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Fig. 4. — The top panel shows high-speed photometry of J1518 
from a representative run, this a portion from 2012 April 16. The 
brightest comparison star is shown in blue, offset by —15%. Av- 
erage point-by-point photometric errors are shown. The bottom 
panel shows a Fourier transform of our entire data set to date, 
more than 32 hr of observations from 2012 March to 2012 July. 
We also display in red the Fourier transform of the residuals after 
pre- whitening by the seven periods listed in the top portion of Ta- 
ble \3\ and mark those periods with green tick marks at the top of 
the panel. We mark the 4(A) significance line as a dashed green 
line. 

TABLE 3 

Frequency solutions for SDSS J151826. 68+065813. 2 



ID 




Period 

(s) 


Frequency 
GuHz) 


Amplitude 
(mma) 


S/N 


h 


2799 


087 ± 0.005 


Family I 

357.2593 ±0.0007 


35 


4 ±0.6 


11.8 


h 


2268 


203 ± 0.004 


440.8777 ± 0.0007 


21 


6 ±0.2 


7.1 


h 


2714 


306 ± 0.003 


368.4183 ± 0.0005 


21 


6 ±0.9 


7.1 


h 


1956 


361 ± 0.003 


511.1532 ±0.0007 


18 


1±0.3 


5.7 


ft 


3848 


201 ± 0.009 


259.8617 ± 0.0006 


15 


7 ±0.3 


5.3 


h 


1335 


318 ±0.003 


748.8855 ±0.0015 


13 


6 ±0.6 


4.5 


h 


2134 


027 ± 0.004 


468.5976 ±0.0008 


11 


2 ±0.4 


4.5 


fi 


2796 


048 ± 0.022 


Family II 

357.6476 ± 0.0028 


11 


1 ±4.8 


13.5 


h 


2209 


996 ± 0.024 


452.4895 ± 0.0049 


19 


9 ±1.4 


6.5 


h 


2802 


769 ±2.6 


356.79 ±0.33 


26 


4 ±3.9 


8.7 


U 


1956 


331 ± 0.008 


511.1609 ±0.0020 


18 


1±0.8 


5.7 


h 


3683 


703 ±7.7 


271.47 ±0.57 


17 


7±3.2 


5.7 


h 


2413 


091 ± 29.5 


414.41 ±5.07 


15 


6 ±4.0 


5.2 


h 


1318 


847 ± 0.004 


758.2382 ±0.0021 


12 


9 ±0.4 


4.4 



mass of the unseen companion is Mi > 0.61 M©, mak- 
ing it most likely another WD. As with J1112, the nature 
of the companion has no direct bearing on pulsations in 
the primary, but it may partially contaminate the spec- 
tral fits. Unlike J1112, no metal lines are detected in the 
spectrum of J1518. 

3.2. Photometric Observations 

We obtained photometric observations of J1518 from 
2012 March to 2012 July for more than 32.2 hr of cover- 
age; a full journal of observations can be found in Table[TJ 
The data were obtained and reduced in an identical man- 
ner to those of J1112, described in Section [2~2l We di- 
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vided the sky-subtracted light curves by the two brightest 
comparison stars in the field, SDSS J151824.11+065723.2 
(g = 14.0 mag) and SDSS J151828.78+065928.9 (g = 
14.3 mag), to allow for fluctuations in seeing and cloud 
cover. 

The top panel of Figure |4] shows a portion of a typ- 
ical light curve for J1518, obtained in 2012 April, and 
includes the brightest comparison star in the field over 
the same period. High-amplitude, multi-periodic, highly 
non-sinusoidal variability is evident. The bottom panel 
of this figure shows an FT utilizing all 16, 522 light curve 
points collected thus far. 

Given the relatively sparse coverage over four months 
(< 1.2% duty cycle), the spectral window for our obser- 
vations of J 15 18 is quite messy, and aliasing from gaps 
in the data makes identifying the underlying periods of 
variability especially difficult. As such, we have provided 
two families of solutions in Table [3l 

The highest peak in the FT occurs at ^357 /iHz. 
However, after pre-whitening by this frequency, there 
is considerable ambiguity as to the next highest peak. 
There are two peaks with nearly the same amplitude, one 
slightly higher at 452.4895 /iHz than the other split apart 
by the daily alias at 440.8777 /iHz. Our two families of 
solutions are based on which of these aliases we pick. 
We note that the top family of solutions holds up more 
robustly to a Monte Carlo simulation of the associated 
frequency and amplitude uncertainties, but the bottom 
family does a marginally better job of reproducing the 
data (the residuals after fitting this solution have a 3.27% 
total amplitude, while the residuals after fitting the top 
family have 3.32% total amplitude). In both cases, the 
fits do a poor job of reproducing the non-linear peakiness 
manifest as the sharp rises and falls in the light curve. 
However, the fits do a decent job of predicting when these 
features occur. 

The seven significant periodicities we identify in J1518 
are by no means exhaustive; we note the significant resid- 
ual power left over after pre-whitening by these periods 
shown as the red FT at the bottom of Figure |U In the 
range between 200 - 800 /Mz (1250 - 5000 s), the aver- 
age amplitude of the original FT is 6.4 mma, and after 
pre-whitening it remains above 3.4 mma. A coordinated 
campaign of near-continuous observations of this rela- 
tively bright WD would greatly help resolve the periods 
present in this pulsating WD. 

4. DISCUSSION 

4.1. Acoustic (p-mode) Pulsations in WDs 

As with the C/O-core DAVs that have been known 
for more than 40 years, the dominant optical variability 
in these ELM WDs is consistent with surface tempera- 
ture variations caused by non-radial g-mode pulsations 
driven to observability by a hydrogen partial ionization 
zone. However, in Jl 112, we also see the first evidence for 
short-period p-mode pulsations in a WD. These p-modcs 
offer a tantalizing opportunity to probe the interior of an 
ELM WD in a complimentary way to the g-modes also 
present. 

Pulsation calculations have long shown that WDs of 
all masses should be unstable to p-mode pulsations 
fg-g. iSaio et~aT1ll983Hs"tarrfield et al.ll983tlHansen et al.l 
1985). However, despite exhaustive searches (e.g. 



IRobinsonlH98l IKawaler et aHfl99i ISilvotti et al1l20H , 
no such high-frequency modes have ever been observed 
in a WD. Discovery of p-mode pulsations in ELM WDs 
could suggest that these pulsation models are indeed 
valid, but that the amplitudes of these oscillations are 
simply too small to detect with significance in the more 
massive DAVs. 

In the context of hot B subdwarf stars (sdBs), it 
is not entirely surprising to find ELM WDs with ob- 
servable p-mode pulsations. Variable hot B subdwarf 
stars (sdBVs), especially of the EC 14026 or V361 Hya 
class, are observed to vary with p-mode pu lsations (e.g. 
Kilkenny et~aL1 Il997t ICharpinet et all 11997ft . These sd- 
BVs are qualitatively quite similar to ELM WDs, with 
He-cores and similar surface gravities, 5.2 < log g < 6.1, 
although they are typically much hotter, with T e g- > 
28,000 K, and more massive, with canonical masses 
^0.45 Mq. Their p ulsation periods are most often in the 
range 120 - 300 s (Kil kenny! [2007ft . quite similar to the 
variability in J1112. Many sdBVs are hybrid pulsators, 
showing both p- and g-mode pulsations. 

Still, we must be careful to rule out all other possi- 
bilities before declaring the short-period instabilities in 
J1112 as b qna-fide p-modes. Th e low-order p-modes cal- 
culated by ICorsico et al.l ()2012ft range from 109 — 7.5 s 
for their 1 < k < 29 models of a 0.17 M He-core WD, 
slightly shorter than the 134.3 s periodicity in J1112. 
However, since that work finds that <?-modes with periods 
shorter than 1100 s are stable to pulsations, p-mode pul- 
sations remain the most likely explanation for the short- 
period variability in J1112. (It also remains possible that 
the 134.3 s periodicity is a radial fundamental mode for 
the star.) 

4.2. Asteroseismic Applications 

While we have not yet matched the observed periods 
in these first three pulsating ELM WDs to full asteroseis- 
mic models, we can begin to assess some of their observed 
properties. We are extremely interested in finding pulsa- 
tion periods that are consecutive radial overtones, since 
the asymptotic period spacings of g-modes are a direct 
probe of the overall mass of the star. We expect that 
low-mass pulsating WDs will have relatively high mean 
period spacings. Fo r the relevant He-core WD models, 
ICorsico et all (|2012ft find that 0.17 M Q WDs have dipole 
forward period spacings of roughly 104 s, compared to 
period spacings of 64 s fora 0.45 Mq He-core WD. Sim- 
ilarly, iSteinfadt et~aTl (|2010ft find an 89 s mean period 
spacing for a 0.17 M Q WD. 

The rich and relatively simple pulsation spectrum of 
Jl 1 12 offers an excellent opportunity to search for con- 
secutive radial overtones. The two shortest-period g- 
mode pulsations in this star, at 1884.6 s and 1792.9 s, 
differ by 91.7 s. If these are i — 1 modes of consecutive 
radial order, this would be excellent direct confirmation 
that this is indeed a roughly 0.17 M Q He-core WD. How- 
ever, a full asteroseismic investigation is warranted to 
confirm that these are indeed consecutive modes. That 
hinges on constructing enough He-core WD models with 
sufficiently different hydrogen layer masses. 

Additionally, J1112 (relatively sinusoidal) and J1518 
(extremely nonlinear) offer strongly contrasting pulse 
shapes. It has been suggested that the large nonlinear 
distortions in the emergent flux are due to the changing 
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TABLE 4 

Properties of the Three Known Pulsating ELM WDs 



Property 



Value 



Property 



Value 



Fig. 5. — The extended ZZ Ceti instability strip. The 56 previously known C/O-core DAVs characterized by Gianninas ct al. (2011) are 
included as purple dots, and our new pulsating ELM WDs are marked in burnt orange. We denote an extrapolated theoretical blue edge 
for the low-mass ZZ instability strip; this dotted blue line is described in the text and u ses the criterion 2tttq = 100 s, with the convective 
prescription ML2/«=1.5. We include the theoretical blue edge for low-mas s DAV s from|Van Grootcl ct al. ( 2013}) as a dashed-dotted blue 
line. We also mark the empirical blue- and red-edges from Gianninas ct al. (2011) as dashed lines. Objects not observed to vary to at least 
greater than 10 mma (1%) are marked with Xs. W e inc lude three new WDs n ot observed to vary, listed in Table [5] the others from past 
work, in gray, were detailed in Hermes ct al. (2012) and Stcinfadt et al. (2012). 

depth of the star's convection z one during a pulsa tion 
cycle (IBrickhilll [l99l IWul l200ll; I Montgomery! [20051) . I f 
this is the case, the method of lMontgomerv et al.l f|2010h 
could be used on J1518 to infer the thermal response 
timescale and thereby the average depth of the convec- 
tion zone in this highly non-sinusoidal DAV. Distortions 
in the light curve may also be the result of other pro- 
cesses, such as the nonlinear resp onse of the emergent 
flux t o temperature perturbations (jFontaine fc Brassard! 
120081) . 

Finally, we can put all three pulsating ELM WDs 
into context by calculat ing their weighte d mean periods 
(WMPs), as defined by iMukadam et al.l (|2006[ ). In the 
classical C/O-core ZZ Ceti stars, cooler WDs typically 
have longer WMPs. This generally holds for the pulsat- 
ing ELM WDs, as well: The coolest of the three, J1840, 
has a WMP -3722 s. J1112 and J1518 have WMPs 
—2288 s and —2404 s, respectively. By comparison, no 
classical C/O-core DAV known before the discovery of 
pulsating ELM WDs had a WMP > 1200 s (or T off < 
10,000 K). 

These pulsating ELM WDs lie in a new region of the 
classical ZZ Ceti instability strip, extending the strip to 
much cooler and lower surface gravity WDs. We have 
used our new temperature and surface gravity determina- 
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tions described in Section ["2. 1.1 1 to characterize our three 
new ELM WDs. This means we have also updated the at- 
mospheric parameters for J1840, the first pulsating ELM 
WD discovered. We also outline th e new period determ i- 
nations from J1840, as detailed in iCorsico et al.l (|2012D . 
These values can be found in Table [U 

We continue to search for new pulsating ELM WDs, 
but there is also utility in constraining the regions 
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TABLE 5 

Observed Low-Mass DAV Candidates and Null Results 



Object 


go-SDSS 


T e g 


log g 


Reference 


Det. Limit 




(mag) 


(K) 


(cm s -1 ) 




% 


SDSS J012549.37+461920.1 


15.8 


9050 ± 580 


7.18 + 0.12 


Brown et al. (2012) 


0.3 


SDSS J144342. 74+150938.6 


18.6 


8810 ± 320 


6.32 + 0.07 


Brown et al. (2(5T2) 


0.1 


SDSS J221928. 48+120418.6 


17.7 


9700 ± 160 


5.12 + 0.08 


Brown et al. (2012) 


0.3 



of parameter space where WDs are not observed to 
vary. In addi t ion to the targets shown not to vary by 
iHermes et al.l (|2012t) . we have put limits on the lack 
of variability in three new ELM WDs, detailed in Ta- 
ble p Our limits on SDSS J144342.74+150938.6 are 
particularly stringent since we observed its field for 
more than 7.2 hr; one of the comparison stars we used, 
SDSS J144347.31+150841.9, turned out to be variable, 
most likely a <5-Scuti star, dominated by a 3.46-hr pul- 
sation period. We construct an updated empirical insta- 
bility strip using the three new pulsators and more than 
a dozen low-mass WDs not observed to vary, shown in 
Figure [5] 

The addition of two new pulsating ELM WDs aids 
in constraining t he theoret i cal low -mas s WD instability 
strip Following IBrickhilll (|1991D and IGoldreich fc Wul 
(1999), we use the criterion that P max ~ 27ttc for the 
longest period mode that is excited, where P max is the 
mode period and the timescale rc describes the heat ca- 
pacity of the convection zone as a function of the local 
photospheric flux, which we compute from a grid of mod- 
el^. Using the traditional value of P max = 100s, we 
find that the convection parameters ML2/a = 1.5 pro- 
vide a good match to the blue edge for the normal-mass 
DAVs, but that extrapolating this relation to the pul- 
sating ELM WDs yields a blue edge which is somewhat 
hotter than the observati ons (see the dotted b lue line in 
Figure EJ). The models of ICorsico et~aTI (|2012[ ) find that 
only g-mode periods greater than ^1000 s are excited, 
but using this criterion only moves the blue edge a small 
amount closer to the location of the observed pulsators 
(about 130 K for a given surface gravity). It is possi- 
ble there are additional mechanisms at work that move 
the blue edge to cooler temperatures. Using a less effi- 
cient prescription for convection, such as ML2/a = 1.0, 
will shift the blue edge to cooler temperatures, as seen in 
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iVan Grootel et all 120131 whose blue edge we include in 
Figure [3 We are also still dominated by small number 
statistics, and perhaps we have simply not yet found any 
hotter pulsating low-mass WDs. It remains difficult to 
predict a theoretical red edge for even the classical C/O- 
core DAVs, so we have not extended it to low-mass WDs, 
an d instead provide a port ion of the red-edge established 
bv lGianninas et al.l (|201lD . 

5. CONCLUSIONS 

We have discovered the second and third pulsating 
ELM WDs as part of a search for such objects at the 
McDonald Observatory. Spectroscopic fits to both ob- 
jects indicate each has a mass < 0.25 M Q , and thus both 
likely harbor a helium core. These variable WDs are 
most likely members of a low-mass extension of the clas- 
sic ZZ Ceti instability strip. 

We will continue our efforts to populate the empir- 
ical low-mass WD instability strip, and look forward 
to the exciting prospect of performing asteroseismology 
by matching the observed pulsation periods in low-mass 
WDs to theoretical He-core WD models. 
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